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i-C ■ Abstract 

^ ' We study the photon statistics in the superpositions of coherent states la) and \a* 

named "Schrodinger real and imaginary cat states". The oscillatory character of the 
photon distribution function (PDF) emerging due to the quantum interference between 
the two components is shown, and the phenomenon of the quadrature squeezing is 
observed for the moderate values of \a\ ~ 1. Despite the quantity (An 2 )/(n) tends 
to the unit value (like in the Poissonian PDF) at \a\ 3> 1, the photon statistics is 
essentially non-Poissonian for all values of \a\. The factorial moments and cumulants 
of the PDF are calculated, and the oscillations of their ratio are demonstrated. 
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1 Introduction 

The coherent states \a) introduced by Glauber |[| with respect to problems of quantum 
optics are eigenstates of the photon annihilation operator a (a\a) = cx\a)), their eigenvalues 
covering the entire complex plane. Since the Heisenberg uncertainty relation is minimized 
in coherent states, the latter are interpreted as the "classical" ones. The photon distribution 
function (PDF) in the coherent state is the usual Poisson distribution, which is a smooth 
function with a single maximum near the mean photon number. 

In nonclassical Gaussian states of light the PDFs possess a number of specific properties 
related to the effect of quadrature squeezing || ||. The PDF for one- mode Gaussian states 
of electromagnetic field was discussed in f5|. In the single- mode correlated state of light 
was introduced, which minimizes the Schrodinger - Robertson uncertainity relation [^|, §]. 
The oscillatory character of PDF for squeezed and correlated Gaussian light was discussed 



in [0, 10 1 . Different nonclassical states of light might play an essential role in nonlinear 



processes 11 . 



By using different finite or infinite sums of coherent states for a one-mode quantum 
harmonic oscillator various non-Gaussian nonclassical states of light may be constructed. 
Recently, there have been much interest in the properties and in the generation of so-called 



"Schrodinger cats" |L2| in the context of quantum optics [^3|| -|p!5|. These quantum su- 
perpositions of the finite or infinite number of coherent states have various nonclassical 
characteristics emerging due to the quantum interference between summands. 



The first example of such kind was constructed in [16[ , where the concept of the even and 



odd coherent states |a±) = iV±(|ot|)(|a) ± | — a)) was introduced. These states satisfy the 
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equation a 2 \a±) = a 2 \a±). Moreover, they manifest such remarkable nonclassical properties 
as oscillations in photon number distribution, sub-Poissonian statistics, and the quadrature 
squeezing. The methods of generating and detecting the even and odd states in quantum 



optics were discussed, e.g., in |L7|, [Lq| . 

The aim of the present article is to study the photon statistics in another representative 
of nonclassical states: the "real" (RCS) and "imaginary" (ICS) coherent states (or "real" 
and "imaginary cats"). 

\^) = N ± (a,a*)(\a)±\a*)), (1) 
the normalization constant being given by (a = |a|e l¥> ) 

N±(a,a*) = ^(l±Ree a2 - |Q|2 r 1/2 = (l ± e- 2|a|2sin2 ^cos(|«| 2 sin2^)) _1/2 . (2) 



y/2 v 1 V2 

Similar superpositions, \a) + \ — a*), were considered from the group-theoretical point of 
view in |19| , where they were called "charged Schrodinger cats" . All the states of this family 
belong to the class of superposition states discussed in |T3| and [ED]. The interference between 



| a) and \a*) leads again to the non-Poissonian (modulated Poissonian) PDF. Evidently, the 
statistical properties of real and imaginary cat states, as well as the degree of squeezing, 
depend essentially on both the modulus |a| and the phase difference 2ip between a and a*. 
It is interesting, however, that the correlation coefficient between the quadrature components 
identically equals zero. 

We obtain the generating function G(n) of factorial moments F(n), and use it to find 
the analytical expressions for the photon means, dispersions, and Mandel's Q-parameter. 
Besides, we calculate the cumulants K{n) of the PDF and demonstrate the oscillatory be- 
haviour (essentially nonclassical) of the ratio H(n) = F(n)/K(n). 



The parameter H(n) was shown in [21 to play an essential role in analysing the particle 



multiplicity distributions in high energy physics, since in experiments it demonstrates an 
oscillatory behaviour, which is sensitive to the detailes of the state of particles created in 



high energy collisions p^| . For the squeezed and correlated states of the electromagnetic 



field this parameter was discussed in [23|. In the present article we concentrate on studying 



the behaviour of H(n) for the Schrodinger real and imaginary cat states. 

2 Means and variances of quadratures 

Using the well known expansion 



oo n 

l«> = e-^ |2 £^|n>, (3) 

n=0 



it is easy to derive the corresponding expansions for RCS and ICS: 



oo I „ \n 

1 1 i2 a 



1ST) = 2iV + e-^l J2 V=f C0S MN> W 
n=o Vnl 



|ST) = 2iN„e->\ ^7=sin(n^)|n), (5) 
n=i vnl 

where < (p < 2ir, and |n) is the vector of the Fock space. The average values of the ladder 
operators a, a' in real and imaginary cat states are invariant with respect to the Hermitian 
conjugation: 

(a)± = (a f )± = 2|iV±| 2 (Re a ± Re («e Q2 " |a|2 )) . (6) 
Consequently, the quadrature operators p, q (we assume H = 1) 

a. — a? A a + a) 
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possess the following mean values: 

(q)± = V2(a ± ), (p)± = 0. (8) 
For the elements of the quadrature covariance matrix we get 

<r(J) = -- 8|iV±| 4 (Ima) 2 (|e a2 -l Q l 2 | 2 ± Ree^H 2 ) , (9) 
°$ = \ + 4|iV±| 2 (Ima) 2 , (10) 

^? = ^<P9 + 5P)± - (P)±<?)± = 0. (11) 

The last equality shows that both real and imaginary Schrodinger cat states have no cor- 
relation between the quadratures. Being rewritten in terms of |a| and ip, the quadrature 
variances read 

a W = X - - 8 (jaiV± 2 | sinv9e- |Q|2sin ^) 2 (e~ 2|a|2sin2 ^ ± cos(|a| 2 sin(2v9))) , (12) 

<rg) = i + 4(|aiV±|sinv9) 2 . (13) 

The variance of the quadrature p always exceeds |. As to the quadrature g, it is squeezed 
provided the inequality 

exp(— 2|a| 2 sin 2 </?) ± cos(|o;| 2 sin2</?) > 

holds. Assuming \a\ 2 to be small, one can expand the summands in powers of \a\ 2 . For the 
real cat state the limit of the inequality is 2 + O (\a\ 2 ) > 0, for any phase (p. The analogous 
limit for the imaginary cat state is less than zero for arbitrary phase. Thus, there is no 
quadrature squeezing in the ICS when \a\ 2 tends to zero. On the contrary, the squeezing in 
RCS exists for any ip up to \a\ 2 = ir/2. 
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3 Photon distribution function, moments and cumu- 
lants 

Using expansions (Eh and (H) one gets the photon distribution functions 



2 a 



2h 



P±(n) = |(n|c^>| 2 = 2|i\r ± (|a|^)| 2 e-l a l i-L-(l ± cos(2r^)). (14) 



TV. 



It is obvious that P±(n), being a modulated Poissonian distribution, is the sum of three 
Poissonian distributions with complex means, which became the usual Poissonian distribu- 
tion with (n) = a 2 at real a. Oscillations of PDF are demonstrated in Figs. 1 and 2 both 
for RCS (a) and ICS (b) in the representative cases of \a\ 2 = 1.5 (Fig. 1) and \a\ 2 = 9.0 
(Fig. 2), ¥? = 0.3vr. 

The simplest way to derive the statistical characteristics of the field (in particular, (n) 
and (An 2 ) = (n 2 ) — (n) 2 ) is to differentiate the generating function 

oo 

G±{\) = P±(n)X n = 2|iV ± | 2 e- |a|2 (e |o|2A ± Ree a2A ) (15) 

n=0 

with respect to the auxiliary real parameter A. Then the factorial moments read 
d n G 



F(n) ee 



n 



= 2|AT ± (|«| 2 ,^)| 2 |«| 2n (l ± e- 2|a|2sin2 ^ cos(|«| 2 sin2^ + 2mp)) . (16) 

A=l 

The factorial moments up to the 10-th order are plotted in Fig. 3 for RCS (a) at \a\ 2 = 1.5 
and (p = 0.37T. The plot of F(n) for ICS at the same a seems similar to Fig. 3. 

It is clear that in the case of \a\ 2 sin 2 cp 1 expression fllfi]) tends to the classical (Poisso- 
nian) limit \a\ 2n . Nevertheless, the PDFs in both RCS and ICS dramatically differ from the 
Poissonian distribution for all values of \a\. This means, in particular, that the knowledge 
of a few of the lowest factorial moments is not sufficient to decide whether the statistics is 



Poissonian or not. In the limit \a\ ^> 1 the difference becomes clear only for the moments 
having the orders not less than n! ks \a\ 2 sin 2 ip/ In \a\. In the case of \a\ 2 = 9.0 and p> = 0.37T 
we have n' ~ 6. 

The mean photon number equals 

dG 



(n)-_ 



A=l 



2|a| 2 |iV ± (|a| 2 ^)| 2 (l±e- 2|a|2sinV cos(|a| 2 sin2v9 + 2v9)) , (17) 



d\ 

whereas Mandel's Q-parameter 

(An 2 ) ± - (n)± 



Q 



± 



(n) z 
is equal to 



Q± = _ \2aN M _ e _ 2H2 sin2 v gin2 , g2 e _ 2|a|2 sin2 „ ± 2 g . n ^ + x 

(n)± v ' ' 

However, in the case under study this parameter is not suitable for the analysis of the 
photon statistics, since for large values of |asin</?| it is close to zero, although statistics 
remains essentially non-Poissonian. Actually the equality Q ~ is only necessary, but not 
sufficient condition of the Poissonian statistics. 

More informative characteristics of the photon statistics are the cumulants 

d n \nG 
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Kin) 



A=l 



To find the relation between the moments and cumulants note that 

d dlnG 

d\ G - G ^>T- 

Thus, 



(19) 



d n _ d n - 1 ( d\nG \ _ y^V k ) d k+1 In Gd n ^- k G 



Putting A = 1 in ( P0| ) one gets the recursive relation 

^ = t { J n m Z k y K ^ n -^ < 21) 

followed by another one, which is useful for the numerical calculation of the cumulant of the 
n-th order through the cumulants of preceding orders and moments up to the n-th order, 

In particular, 

Kt = F 1 , 

K 2 = F 2 -Ff, 

K 3 = F 3 -3F 2 F l +2F?, 

K 4 = F 4 - 4F 3 F 1 - 3F 2 + 12F 2 F 2 - 

K 5 = F 5 - WF 2 F 3 - 5FxF 4 + 20F 2 F 3 + 30F^F 1 - 60FfF 2 + 24F X 5 , 

K 6 = F 6 - 6FxF 5 - 15F 2 F 4 + 120F 1 F 2 F 3 + 30F 2 F 4 - 120F X 3 F 3 
-270F|iJ? + 360i^F a - 10F 2 + 30F 2 3 - 120-Ff. 

The cumulants of the photon number up to n = 10 corresponding to |5^T)-state with 
| ck | 2 = 1.5, <f = 0.37T are plotted in Fig. 4 (a) to demonstrate the fast (factorial-like) diver- 
gence of a series of K(n). Since the cat states with \a\ 2 ~ 1 are inherently nonclassical, 
the cumulants differ from zero even for small values of n. Their oscillating behavior is 
demonstrated in Fig. 4 (b). 
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The ratio H(n) = F(n)/K(n) may be used instead of Q-parameter to evaluate how 
"nonclassical" is a quantum state |23 |. Its oscillatory behavior in RCS and ICS with \a\ 2 = 
1.5 and \a\ 2 = 9.0, ip = 0.3n is seen from Fig. 5. The strong decrease at higher n is also 
seen: the factorial-like divergence of K(n) surpasses the quasiexponential growth of F(n). 

The high amplitude oscillations of H(n) = F(n)/K(n) for cat states with | ct | 2 = 9.0 
admit a simple explanation : while F(n) equals approximately the classical value |a| 2n , 
the cumulants K(n) pa ±0. Nevertheless, for large n the ratio F(n)/K(n) damps even 
for | ct | 2 = 9.0. Note that the plot of H(n) for the imaginary cat state with | ct | 2 = 9.0 is 
symmetrical to the analogous plot for the real cat state with respect to n-axis. 



4 Conclusion 

We have shown that both real and imaginary Schrodinger cat states manifest quadrature 
squeezing and non-Poissonian statistics. The factorial moments of the photon distribution 
function have proven to be monotonically growing, while the oscillatory factorial divergence 
to infinity is observed for the cumulants K(n) when n>l. The ratio F(n)/K(n) has been 
found to be strongly oscillating for small n and damping when n grows. These properties 
confirm the nonclassical nature of the Schrodinger real and imaginary cat states. 

If one considers the time evolution of quantum states in the case of a harmonic oscillator 
with unit frequency, then RCS and ICS will be transformed as 

|o±Tt) = N±(a,a*) (\ae~ u ) ± |a*e^*)) 

(with the same complex exponentials in both terms). Thus in the time- dependent case the 
complex parameters of the coherent states determining the superpositions do not preserve 
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the property to be mutually complex conjugated. However, the phase difference 2ip between 
these parameters is time invariant. Since the PDF is also time invariant, it is convenient 
and natural to choose the initial moment in such a way that the centers of the two coherent 
states wave packets would be located symmetrically with respect to the real axis. Just this 
natural choice has been made in the paper. 

The multimode generalization of the real and imaginary cat states may be constructed 



analogously to the multimode generalization of the even and odd coherent states done in [23 . 
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Figure captions 

Fig. 1. Photon distribution functions P+(n) for |5qT) (a) and P-(n) for \aZ) (b) with \a\ 2 = 
1.5, ip = 0.3n. 

Fig. 2. P + (n) (a) and P_(n) (b) for \a\ 2 = 9.0, = 0.3vr. 

Fig. 3. Factorial moments F(n) of the photon number for |<xjT) (a) and \ocl) (b) with 
\a\ 2 = 1.5, up = 0.3tt. 

Fig. 4. Cumulants if (n) of the photon number for |SZjT)-state up to the 10-th order (a) and 

up to the 5-th order (b) with \a\ 2 = 1.5, (p = 0.3n. 
Fig. 5. Ratio F(n)/K(n) for \a^) (a,c) and \aZ) (b,d) with \a\ 2 = 1.5 (a,b), \a\ 2 = 9.0 (c,d) 

and ip = 0.3n. 
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